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ABSTRACT 

Using new and published data, we construct a sample of 160 brightest cluster galaxies 
(BCGs) spanning the rcdshift interval 0.03 < z < 1.63. We use this sample, which 
covers 70% of the history of the universe, to measure the growth in the stellar mass of 
BCGs after correcting for the correlation between the stellar mass of the BCG and the 
mass of the cluster in which it lives. We find that the stellar mass of BCGs increase 
by a factor of 1.8 ± 0.3 between z — 0.9 and z = 0.2. Compared to earlier works, our 
result is closer to the predictions of semi-analytic models. However, BCGs at z = 0.9, 
relative to BCGs at z = 0.2, are still a factor of 1.5 more massive than the predictions 
of these models. Star formation rates in BCGs at z ~ 1 are generally too low to result 
in significant amounts of mass. Instead, it is likely that most of the mass build up 
occurs through mainly dry mergers in which perhaps half of the mass is lost to the 
intra-cluster medium of the cluster. 

Key words: galaxies: clusters: general - galaxies: evolution - galaxies: high-redshift 
- cosmology: observations 



1 INTRODUCTION 

Brightest Cluster Galaxies (BCGs) are amongst the largest, 
most luminous and most massive galaxies in the universe at 
the present epoch. Located in the cores of rich galaxy clus- 
ters, BCGs are easy to identify, both observationally and 
in simulations. They can also be observed at a time when 
the universe was less than a third of its current age. They 
therefore provide an attractive target for testing our under- 



* E-mail: clidman@aao.gov.au 
© 2002 RAS 



standing of the processes that drive galaxy evolution, albeit 
in the most massive galaxies of the universe. 

In the hierarchical scenario for the formation of struc- 
ture in our universe, galaxies start off as small fluctuations 
in the density of matter and build up their stellar mass 
over time by converting material accreted from their sur- 
roundings into stars and by merging with other galaxies (see 
Baugh 2006, for a review). In semi-analytic models that use 
the hierarchical scenario as their foundation, the stellar mass 
of a BCG increases significantly with time. For example, be- 
tween redshift z = 1.0 (corresponding to a look-back time of 
6.7 Gyr) to z = 0, the semi-analytic model described in De 
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Lucia & Blaizot (2007) predicts that BCGs increase their 
stellar mass by a factor of four (De Lucia & Blaizot 2007). 

In contrast to this prediction, observations appear to 
suggest that there is little growth in the stellar mass of 
BCGs, although apparently conflicting results have been 
reported. Using a sample of optically selected clusters, 
Aragon-Salamanca, Baugh, & Kauffmann (1998) found that 
the stellar mass of BCGs grew by a factor of 4 between 2 = 1 
and today. Burke, Collins, & Mann (2000), on the other 
hand, using a sample of X-ray selected clusters over a sim- 
ilar redshift range, find substantially less growth. 1 Burke, 
Collins, & Mann (2000) conclude that sample selection can 
explain part of the difference between their results and those 
in Aragon-Salamanca, Baugh, & Kauffmann (1998), a con- 
clusion that was supported by Nelson et al. (2002). In an 
independent study, using an optically selected sample of 21 
high-redshift clusters, Whiley et al. (2008) find little change 
in the stellar mass of BCGs since z ~ 1. 

At higher redshifts, the discrepancy between the mod- 
els and the observations is larger. Collins et al. (2009) and 
Stott et al. (2010), using a sample of 20 mostly X-ray se- 
lected clusters and a sample of nearby clusters from Stott et 
al. (2008), find that there is little growth between z ~ 1.4 
and now. At z ~ 1.4, the semi-analyic model of De Lucia & 
Blaizot (2007) predicts that BCGs should be a factor of six 
less massive. Therefore, there appears to be clear disagree- 
ment between the models and the observations. 

In this paper, we expand upon the work that has been 
done so far in three ways. First, we increase the number of 
BCGs beyond z = 0.8 for which accurate near-IR photome- 
try is available. Second, we extend the redshift baseline by 
including the BCGs in two recently discovered clusters at 
z ~ 1.6. Third, we use our large sample to account for the 
correlation between the stellar mass of the BCG and the 
mass of the cluster in which it lives. 

We start the paper with a description of our new sam- 
ple of BCGs in Section 2, followed, in Section 3, with a 
description of the near-IR imaging data that we use in 
later sections. In Sections 4 and 5, we derive the magni- 
tudes and colours of the BCGs in our sample and compare 
them to predictions made by simple and composite stellar 
population models. Following Stott et al. (2010), we use 
this comparison to estimate stellar masses. In Section 6, 
we discuss our results, comparing them to the predictions 
made by semi-analytic models and examining how robust 
they are to our methods. In the final section, we summarise 
our main results. Throughout the paper, all magnitudes and 
colours are measured in the observer frame and are placed 
on the 2MASS photometric system. Vega magnitudes are 
used throughout the paper. We also assume a flat ACDM 
cosmology with Aa = 0.73 and Hq — 70 kms -1 Mpc -1 . 



2 A NEW SAMPLE OF DISTANT BCGS 

We use clusters from the SpARCS 2 survey (Muzzin et al. 
2009; Wilson et al. 2009; Demarco et al. 2010; Muzzin et al. 
2012) to assemble a sample of 12 BCGs spanning the red- 
shift interval 0.85 < z < 1.63. The coordinates and redshifts 
of the clusters are listed in Table 1. Ten of the twelve clus- 
ters were observed in the GCLASS 3 survey, which used the 
Gemini Multi-object Spectrographs on Gemini North and 
Gemini South Telescopes to obtain between 20 and 80 spec- 
troscopically confirmed members per cluster (Muzzin et al. 
2012). The other two clusters, which are the most distant 
clusters in our sample, are more recent discoveries. Both 
clusters are spectroscopically confirmed, with a dozen spec- 
troscopic redshifts per cluster (Muzzin et al. in preparation; 
Wilson et al. in preparation). 

All twelve clusters were discovered by searching for over- 
densities in the number of red galaxies using a combina- 
tion of images taken with IRAC on the Spitzer Space Tele- 
scope with z-band images taken with either MegaCam on 
the Canada-France-Hawaii Telescope (CFHT) or MOSAIC 
II on the Cerro Tololo Blanco Telescope. The ten GCLASS 
clusters were found using the z-[3.6] colour, whereas the two 
more distant clusters were found using the [3.6]-[4.5] colour 
together with the requirement of a red z-[3.6] colour. Further 
details of how the clusters were discovered can be found in 
(Muzzin et al. 2008, 2009; Wilson et al. 2009) and (Muzzin 
et al., in preparation). 



3 OBSERVATIONS 

We used three near-IR imaging cameras to observe 12 clus- 
ters. Six of the clusters were observed with the Wide-field 
InfraRed Camera (WIRCam) on the Canada-France-Hawaii 
Telescope (CFHT) on Mauna Kea. Another three clusters 
were imaged with the Infrared Side Port Imager (ISPI) on 
the Blanco Telescope at the Cerro Tololo Inter- American 
Observatory (CTIO). Finally, the three most distant clusters 
were imaged with the High Acuity Wide field K-band Imager 
(HAWK-I) on Yepun (VLT-UT4) at the ESO Cerro Paranal 
Observatory. The fields-of-view of the imagers, and their 
plate scales are noted in Table 2. Details of the observations, 
including exposure times, are listed in Table 1. With the ex- 
ception of the two most distant clusters (SpARCS J033056- 
284300 and SpARCS J022426-032331), all clusters were im- 
aged in J and Ks. At the time SpARCS J033056-284300 
and SpARCS J022426-032331 were observed with HAWK-I, 
neither cluster had been spectroscopically confirmed. The 
Y and Ks filter pair were chosen over J and Ks, since the 
former pair almost straddle the 4000 A break, thereby in- 
creasing the contrast of cluster members over field galaxies 
and easing target selection for spectroscopy. 



1 Both Aragon-Salamanca, Baugh, & Kauffmann (1998) and 
Burke, Collins, & Mann (2000) use an Einstein de-Sitter universe, 
i.e. f2Mi^A = 1,0, with Ho = 50 kms -1 Mpc~ 1 for the cosmol- 
ogy. While their results are not directly comparable to the results 
in later papers, one can compare the results of the two papers. 



2 Spitzer Adaptation of the Red-Sequence Cluster Survey, 
www.faculty.ucr.edu/~gillianw/SpARCS/ 

3 Gemini Cluster Astrophysics Spectroscopic Survey, 
www.faculty.ucr.edu/~gillianw/GCLASS/ 
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Table 1. Observational summary 



Cluster 


Redshift 


R.A. 5 


Decl. 5 


Instrument /Telescope 


Exposure times: J and Ks 






J2000 


J2000 






M 


SpARCS J003442-430752 1 


0.867 


00:34:42.03 


-43:07:53.4 


ISPI/Blanco 


17280 


8800 


SpARCS J003645-441050 1 


0.869 


00:36:44.99 


-44:10:49.8 


ISPI/Blanco 


17280 


7080 


SpARCS J161314+564930 1 - 2 


0.871 


16:13:14.63 


56:49:30.0 


WIRCAM/CFHT 


6240 


6300 


SpARCS J104737+574137 1 


0.956 


10:47:33.43 


57:41:13.4 


WIRC AM/CFHT 


7560 


2400 


SpARCS J021524-034331 1 


1.004 


02:15:23.99 


-03:43:32.2 


ISPI/Blanco 


26640 


11800 


SpARCS J105111+581803 1 


1.035 


10:51:11.22 


58:18:03.3 


WIRCAM/CFHT 


6840 


2700 


SpARCS J161641+554513 1 - 2 


1.156 


16:16:41.32 


55:45:12.4 


WIRCAM/CFHT 


18960 


7000 


SpARCS J163435+402151 1 ' 3 


1.177 


16:34:38.21 


40:20:58.4 


WIRCAM/CFHT 


11640 


6850 


SpARCS J163852+403843 1 ' 3 


1.196 


16:38:51.64 


40:38:42.8 


WIRCAM/CFHT 


11640 


6000 


SpARCS J003550-431224 1 - 4 


1.335 


00:35:49.68 


-43:12:23.8 


HAWK-I/Ycpun 


11040 


12000 


Cluster 


Redshift 


R.A. 5 


Decl. 5 


Instrument/Telescope 


Exposure times: Y and Ks 






J2000 


J2000 






M 


SpARCS J033056-284300 


1.626 


03:30:55.87 


-28:42:59.7 


HAWK-I/Ycpun 


8880 


3040 


SpARCS J022426-032331 


1.633 


02:24:26.32 


-03:23:30.7 


HAWK-I/Yepun 


8640 


5040 



Note 1: Muzzin et al. (2012) 
Note 2: Demarco et al. (2010) 
Note 3: Muzzin et al. (2009) 
Note 4: Wilson et al. (2009) 
Note 5: Coordinates of the BCG 



Table 2. Instrument summary 



Instrument 


Telescope 


Pixel Scale 
["] 


FoV 

['] 


Detector 


WIRC AM 1 


CFHT 


0.304 


20.5 


2x2 Hawaii-2RG mosaic 


ISPI 2 


Blanco 


0.307 


10.3 


Hawaii-2 


HAWK-I 3 


Yepun (VLT-UT4) 


0.1065 


7.5 


2x2 Hawaii-2RG mosaic 



Note 1: (Puget et al. 2004) 

Note 2: (van der Bliek et al. 2004) 

Note 3: (Pirard et al. 2004; Casali et al. 2006) 



3.1 Data Reduction 

The processing of the raw data was done in a standard man- 
ner and largely follows the steps outlined in Lidman et al. 
(2008). Data from each of the cameras were pre-processed 
(dark subtraction, flat-fielding, and sky subtraction) using a 
combination of observatory-developed instrument pipelines 
(for example, the CFHT data were processed with version 
1.0 of the Tiwi pipeline 4 ) and our own scripts using IRAF 5 . 

SCAMP (version 1.6.2) and SWarp (version 2.17.6) 6 
were used to map the sky-subtracted images onto a common 
astrometric reference frame. After accounting for gain vari- 
ations between chips (only relevant for the data that were 
taken with HAWK-I and WIRCAM) and creating individual 
bad pixel maps to account for bad pixels and remnants from 
bright stars observed in previous frames, the images were 
then combined with the imcombine task within IRAF. Each 



4 http://www.cfht.hawaii.edu/Instruments/Imaging/WIRCam/ 

5 IRAF is distributed by the National Optical Astronomy Ob- 
servatories which are operated by the Association of Universities 
for Research in Astronomy, Inc., under the cooperative agreement 
with the National Science Foundation 

6 http://www.astromatic.net/ 



image was weighted with the inverse square of the FWHM 
of the PSF. 

With the exception of the data taken in the Y band, zero 
points were set using stars from the 2MASS point source 
catalogue (Skrutskie et al. 2006). Typically, between 10 to 
40 unsaturated 2MASS stars with 2MASS quality flags of 
'A' or 'B' were selected to measure zero points and their 
uncertainties. 2MASS stars were weighted by the reported 
uncertainties in the 2MASS point source catalogue. The un- 
certainties in the zero points are generally less than 2%, and 
more typically 1%, for both J and Ks. For Y, the zero point 
was set using standard stars that were observed during the 
same night as the clusters. The uncertainty is estimated from 
the night-to-night variation in the zero points and is around 
2%. 



3.2 Data quality 

Overall, the depth and quality of the imaging data varies 
substantially from one image to another. The image quality, 
as measured from bright stars, varies from 0'.'3 in the data 
taken with HAWK-I to l'/5 in the data taken with ISPI. 

The image depth, which we define as the 5 sigma point 
source detection limit, varies from 19.5 for the Ks band im- 



© 2002 RAS, MNRAS 000, 1-?? 



4 C. Lidman et al. 



age of SpARCS J003645-441050 to 25.1 for the Y band im- 
age of SpARCS J022426-032331. In all cases, the BCG is at 
least 2 mag brighter than the detection limit. The median 
signal-to-noise ratio is around 50. Table 3 summarises the 
image quality and image depth. 



4 ANALYSIS 

Identifying the brightest galaxy 7 in each cluster was gen- 
erally straightforward. Images of the BCGs are shown 
in Fig. 1, and their coordinates are listed in Table 1. 
With only two exceptions - SpARCS J105111+581803 and 
SpARCS J163435+402151 - the BCGs are located near 
to the projected centre of the clusters. For both SpARCS 
J105111+581803 and SpARCS J163435+402151, the BCGs 
are ~ 250 kpc from the projected centre of the cluster. The 
projected distances are not excessively large when compared 
to low-redshift clusters (Bildfell et al. 2008; Sehgal et al. 
2012), and both BCGs have redshifts that place them within 
300 km/s of cluster redshift. 

4.1 Photometry 

To estimate total magnitudes of the BCGs in the SpARCS 
clusters, we follow Stott et al. (2010) and use the SExtrac- 
tor 8 MAG_AUT0 magnitude. MAG_AUT0 is a Kron-like magnitude 
(Kron 1980) within an elliptical aperture. For a given object, 
the elongation and orientation of the aperture is determined 
by second order moments of the light distribution. In this 
paper, the size of the aperture is set to the standard value of 
2.5 times the first raw moment. Note that the definition of 
the first raw moment used by SExtractor differs from the one 
used in (Kron 1980). See the SExtractor user's manual and 
Graham & Driver (2005) for further details. Other SExtrac- 
tror parameters are set to their default values. For example, 
the background is determined globally and the minimum 
Kron radius is set to 3.5, the units of which are not pixels 
but in units of the semi-major (or semi-minor) axis. 

For measuring colours, we first match the image quality 
between images using the IRAF psf match task and then 
measure the flux in apertures that have a diameter of 16 kpc. 
We use the same physical diameter for all BCGs. At z ~ 1, 
this projects to ~ 2'.'1 on the sky. The apertures we use are 
double the size of the apertures used in Stott et al. (2010). 
The image quality of the poorest images - the Ks-band image 
of SpARCS J003645-441050, for example - is not sufficiently 
good enough to use apertures this small. 

We investigated how the colours change with the size of 
the apertures. We varied the aperture diameter from 10 kpc 
to 24 kpc. With three exceptions, the colours change by less 
than 2%, which is similar to the statistical uncertainty. The 
exceptions are the BCGs in SpARCS J105111+581803 and 
our two most distant clusters, where we see changes of up 
to 6%. Interestingly, the BCG of SpARCS J105111+581803 
is about 250 kpc from the centre of the cluster, is an [Oil] 
emitter, and has, relative to other BCGs, a blue colour. Our 

7 Throughout this paper, the BCG is defined as the brightest 
cluster member in the observer-frame Ks band. 

8 We used version 2.6.6 of SExtractor 
http://www.astromatic.net/ - in double image mode. 



two most distant BCGs, which were observed in Y and Ks, 
are also [Oil] emitters. The change in colour with aperture 
diameter might indicate that these galaxies have substantial 
colour gradients. The other BCGs do not show any evidence 
for colour gradients over the range of apertures explored. 

Errors in the photometry are dominated by sky noise, 
so they were estimated by examining the distribution of the 
integrated counts in apertures that were randomly placed 
in regions that were free of objects. For colours, the errors 
were estimated for each filter separately and then added in 
quadrature. 

The filter transmission curves of the J and Ks bands 
in ISPI, WIRCam and HAWK-I are similar to one another; 
however, they differ slightly from the filter curves of the re- 
spective filters used in 2MASS. To account for this difference 
we offset the colours by an amount that depends on i) the 
average spectral energy distribution (SED) of the stars used 
to determine the image zero points, ii) the SED of the BCG 
and iii) its redshift. To determine the offset, we assume that 
the average star can be modelled as a K5 dwarf, which has a 
J-Ks colour that is similar to the average colour of stars that 
are used to determine the zero-point, and that the SED of 
the BCG corresponds to the one predicted by model #3 in 
Fig. 2. See Sec. 5.2 for a detailed description of this model. 
Between redshift 0.8 and 1.6, the magnitude of the correc- 
tion is about 0.08 mag. At these redshifts, the dependence of 
the correction on redshift and the assumed spectrum of the 
BCG is slight, with extreme values differing by 0.03 mag. 

The magnitudes and colours of the 12 BCGs in our sam- 
ple are shown in Table 4. 

4.2 External samples 

To our sample of 12 high-redshift BCGs, we add BCGs from 
a number of external samples. At low to intermediate red- 
shifts (z = 0.04 to 2 = 0.83), we add 103 of the 104 BCGs 
from Stott et al. (2008), excluding the BCG of MS1054.5- 
0321. At higher redshifts (z = 0.81 to z = 1.46), we use 
a sample of 20 BCGs from Stott et al. (2010). The BCG 
of MS1054. 5-0321 is common to both samples. Like the 12 
BCGs in our sample, the BCGs from these two external 
samples were observed in J and Ks. 

For the z < 0.15 BCGs in Stott et al. (2008), Stott et al. 
(2008) used photometry from the extended and point source 
catalogues of 2MASS (Skrutskie et al. 2006). For the z > 
0.15 BCGs, Stott et al. (2008) used the SExtractor MAG_BEST 
magnitude. Depending on the level of crowding, MAGJ3EST 
is either a corrected isophotal magnitude MAG_IS0C0R or the 
MAG_AUT0 magnitude (Bertin & Arnouts 1996). 

The BCGs in the high-redshift sample of Stott et al. 
(2010) are hosted by clusters that come from a number of 
sources. Not all of the clusters are X-ray selected; however, 
all are X-ray luminous, with X-ray luminosities exceeding 
10 44 erg/s. The photometry of these BCGs is measured with 
the SExtractor MAG_AUT0 magnitude. The BCGs in the low- 
to-intermediate-redshift sample of Stott et al. (2008) are all 
hosted by clusters that have X-ray luminosities in excess of 
10 44 erg/s (in the 0.1-2.4 KeV band). 

Additional BCG samples have been published in the 
literature. Aragon-Salamanca, Baugh, & Kauffmann (1998) 
published k-corrected K-band magnitudes for BCGs in 25 
clusters up to z = 0.92. Whiley et al. (2008) combined 
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Table 3. Image quality and image depth 



Cluster 


Image quality 


Image depth 1 


Image quality Ima 


5c dcpth x 




["] 


[mag] 


["] 


mag] 






J 


Ks 




SpARCS J003442-430752 


1.25 


21.8 


0.98 


19.9 


SpARCS J003645-441050 


1.13 


21.7 


1.47 


19.5 


SpARCS J161314+564930 


0.77 


22.2 


0.72 


21.1 


SpARCS J104737+574137 


0.69 


22.2 


0.60 


21.2 


SpARCS J021524-034331 


1.07 


21.8 


0.89 


20.3 


SpARCS J105111+581803 


0.66 


22.5 


0.74 


20.6 


SpARCS J161641+554513 


0.70 


22.8 


0.75 


21.2 


SpARCS J163435+402151 


0.65 


22.9 


0.67 


21.2 


SpARCS J163852+403843 


0.61 


23.1 


0.58 


21.5 


SpARCS J003550-431224 


0.35 


24.6 


0.31 


23.1 






Y 


Ks 




SpARCS J033056-284300 


0.45 


24.0 


0.29 


21.9 


SpARCS J022426-032331 


0.34 


25.1 


0.51 


21.5 



lr The image depth is the 5 sigma point-source detection limit measured over an aperture that has a diameter that is 
twice the image quality. 




SpARCS J003442-430752 , 
^ z=0.867 I 




LI 



Sn^RCS 

* 



J021524-034331 
z=1.004 



* 



SpARCS J163852+403843 
z=1.196 



SpARCS J003 645 -44 1050 
z=0.869 




SpARCS J16 13 14+5 64930 
z=0.871 



* 




SpARCS J105 111+581803 
z=1.035 



SpARCS J161641+5545 13 
z=1.156 



SpARCS J104737+574137 
z=0.956 

% 



SpARCS J163435+402151 



- 



1 > * 3 

Mi 




SpARCS J003 5 5 0-43 1224 SpARCS J033056=#4300 SpARCS J02 



31 



z=1.335 



Figure 1. Ks-band cutouts of the 12 BCGs used in this paper. The images are 9" on a side, which corresponds to 70kpc for the nearest 
BCG and 78kpc for the most distant. 
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Table 4. SpARCS BCG photometry 



Cluster Redshift Ks J-Ks 

[mag] [mag] 



SpARCS 


J003442-430752 





867 


16 


516 


(0.039) 


1.863 


(0.033) 


SpARCS 


J003645-441050 





867 


16 


092 


(0.047) 


1.837 


(0.030) 


SpARCS 


J161314+564930 





873 


15 


693 


(0.015) 


1.794 


(0.012) 


SpARCS 


J104737+574137 





956 


17 


140 


(0.031) 


1.889 


(0.029) 


SpARCS 


J021524-034331 


1 


004 


16 


876 


(0.140) 


1.861 


(0.052) 


SpARCS 


J105111+581803 


1 


035 


16 


877 


(0.046) 


1.740 


(0.030) 


SpARCS 


J161641+554513 


1 


156 


17 


017 


(0.031) 


1.729 


(0.022) 


SpARCS 


J163435+402151 


1 


177 


17 


349 


(0.023) 


1.839 


(0.026) 


SpARCS 


J163852+403843 


1 


196 


17 


647 


(0.052) 


1.913 


(0.051) 


SpARCS 


J003550-431224 


1 


340 


17 


524 


(0.014) 


1.981 


(0.009) 


SpARCS 


J033056-284300 


1 


620 


17 


881 


(0.041) 




(...) 


SpARCS 


J022426-032331 


1 


630 


18 


071 


(0.026) 




(...) 



this sample with 2MASS photometry of the low-redshift 
BCG sample of von der Linden et al. (2007) and their own 
photometry of a sample of 21 intermediate-to-high redshift 
(0.39 < 2 < 0.96) BCGs from the ESO Distant Cluster 
Survey. The photometry of all these samples are measured 
in fixed 37kpc diameter apertures, and is converted to the 
rest-frame K-band using k-corrections. In this paper, as in 
the papers of Stott et al. (2008) and Stott et al. (2010), 
we do not apply k-corrections and we measure the flux in 
differently sized apertures. These differences mean that we 
cannot use the photometry from these studies directly with- 
out inverting the k-corrections and applying a correction for 
the different size of the apertures. Without reanalysing the 
data, the latter is difficult to estimate, so we do not add the 
BCGs from these samples to ours. 

In addition to the BCGs in Stott et al. (2010), we add 
BCGs in 15 X-ray luminous clusters from the intermediate 
redshift CNOC1 9 cluster sample (Yee, Ellingson, & Carlberg 
1996). The CNOC1 clusters are from the Einstein Medium 
Sensitivity Survey (Gioia et al. 1990). We use the Ks-band 
photometry from Muzzin et al. (2007a); we note that these 
clusters were not observed in the J band. 

Ks-band imaging data for 14 of the 15 CNOC1 clusters 
were obtained using the Ohio State-NOAO Infrared Imag- 
ing Spectrograph (ONIS) on the Kitt Peak National Obser- 
vatory (KPNO) 2.1 m telescope. ONIS has a pixel scale of 
0.288", which is similar to the pixel scale of the cameras used 
to observe most of the clusters in our SpARCS sample. One 
cluster, MS 0440+02 was obtained using the PISCES cam- 
era on the Steward Observatory 90 inch (2.3 m) telescope. 
PISCES has a pixel scale of 0.495". The image quality in 
the fully reduced images, varies between 0.7 and 1.3" (See 
Muzzin et al. 2007a, for further details). 

We have reanalysed the processed Ks-band images of 
clusters in the CNOC1 sample following the procedure used 
for clusters in our SpARCS sample (see Section 4.1). The 
Ks-band magnitude of these galaxies is reported in Table 5. 

The BCGs from the four samples (SpARCS, CNOC1, 
Stott et al. (2008), and Stott et al. (2010)) are combined 
into a single sample that is then used to make three sub- 
samples covering three broad redshift ranges: a low-redshift 
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Table 5. CNOC1 BCG photometry 



Cluster Redshift Ks 

[mag] 



A2390 





228 


13 


489 


(0 


068) 


MS0440+02 





197 


13 


337 


(0 


054) 


MS0451+02 





201 


13 


938 


(0 


065) 


MS0839+29 





193 


13 


411 


(0 


062) 


MS1006+12 





261 


13 


786 


(o 


074) 


MS1231+15 





235 


13 


891 


(0 


065) 


MS1455+22 





257 


13 


558 


(0 


062) 


MS0016+16 





547 


15 


288 


(0 


078) 


MS0302+16 





425 


15 


008 


(0 


065) 


MS0451-03 





539 


15 


176 


(0 


071) 


MS1008-12 





306 


13 


676 


(o 


076) 


MS1224+20 





326 


14 


409 


(0 


078) 


MS1358+62 





329 


14 


292 


(o 


063) 


MS1512+36 





373 


14 


632 


(0 


082) 


MS1621+26 





427 


14 


977 


(o 


069) 



subsample (0.0 < z ^ 0.3), an intermediate-redshift sub- 
sample (0.3 < z ^ 0.8) and a high redshift subsample 
(0.8 < z < 1.65). The number of BCGs in each of these 
subsamples is listed in Table 7. We use these subsamples 
throughout the rest of this paper. 



5 ESTIMATING THE STELLAR MASS 

Following the methods used in previous works on determin- 
ing the stellar mass of BCGs (Stott et al. 2008; Collins 
et al. 2009; Stott et al. 2010), we use the offset between 
the observed and predicted observer-frame Ks-band mag- 
nitudes to estimate stellar mass. The predicted magnitude 
is estimated from stellar population models that match the 
observer-frame J-Ks colour of the BCGs over the entire red- 
shift range covered by our subsamples, i.e. from z ~ to 
z ~ 1.6. When converting between luminosity and stellar 
mass we assume that the mass-to-light ratios of the BCGs 
are independent of stellar mass. 
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5.1 Modelling the J-Ks colour 

We use the simple stellar population (SSP) models from 
Bruzual & Chariot (2003, hereafter BC03) to model the 
evolution of the spectral energy distributions (SED) with 
cosmic time. There are a number of ingredients that go into 
the models, such as the initial mass function (IMF), the age 
and duration of the star burst, the metallicity of the stars, 
and dust extinction. 

We assume a Chabrier IMF (Chabrier 2003). The differ- 
ence in the resulting J-Ks colour from using a different IMF 
(for example, a Salpeter IMF (Salpeter 1955)) is less than 
0.02 mag over the entire redshift range covered by the obser- 
vations. Similarly, the stellar mass ratio between BCGs at 
low and high-redshift is relatively unaffected by our choice 
of the IMF. The stellar masses themselves, however, change 
significantly. Excluding stellar remnants, the difference is 
about a factor of two for a given Ks-band luminosity. In this 
paper, we do not use the stellar masses directly, just their 
ratios. 

We assume that extinction from dust in negligible. From 
the small amount of scatter in the colour of galaxies on 
the red-sequence, one can infer that dust either reddens all 
galaxies by a small amount or reddens a small number of 
galaxies considerably (Meyers et al. 2012). If the former is 
true, then the amount of reddening affecting the BCGs in 
our sample is unimportant. If the latter is true, then we 
would expect to see significant colour outliers in Fig. 2, 
which we do not see. 

With the IMF set and dust ignored, we consider a se- 
ries of models in which we allow the star formation his- 
tory to vary. We add an extra dimension to these models 
by combining two models with identical star-formation his- 
tories but different metallicities: a solar metallicity model 
and a model that is two-and-a-half times solar. At low red- 
shifts (z ~ 0.03), BCGs have metallicities that are around 
twice solar (Loubser et al. 2009). We add this extra degree 
of freedom because it is not possible to match the colours of 
the BCGs over the entire redshift range - even by varying 
the star-formation history - with the range of metallicities 
available in BC03. We allow the mass ratio of the two com- 
ponents to vary over the full range (i.e. 0.0 to 1.0) in steps of 
0.1. We also allow the e-folding time, r, of the star-formation 
rate to vary between 0.3 and 1.0 Gyr in steps of 0.1 Gy. We 
then choose the model that best fits the data by finding the 
model with the smallest chi-square. With the exception of 
the CNOC1 sample and the two most distant clusters in the 
SpARCS sample, which lack J band data, we use the en- 
tire sample when fitting the models. If an error in the J-Ks 
colour is unavailable, we assume and error of 0.1 mag. 

The best fitting model, model 3 in Fig. 2, has an e- 
folding time of r = 0.9 Gyr and a composition that is 
split 60/40 between solar metallicity and a metallicity that 
is two-and-a-half times higher. This model accurately de- 
scribes the general change in J-Ks colour with redshift al- 
though it does not capture the scatter. In the remainder of 
this paper we use this model to estimate stellar masses. 

In order to demonstrate the sensitivity of the J-Ks 
colour on metallicity and different star formation histories, 
we plot a series of models. The parameters defining the mod- 
els are listed Table 6. In model 1, all the stars form in a 
single burst at z = 2. In this model, all the stars have solar 



metallicity. Model 2 is similar to model 1, except that we 
move the burst to z — 5 and increase the metallicity to two- 
and-a-half times solar. In model 4, we move away from a 
single burst, using instead an exponentially decaying burst 
of star formation with r = 1.03 Gyr that starts forming stars 
at z = 10.6. This model mimics the star formation history 
of BCGs in the hierarchical models of De Lucia & Blaizot 
(2007), in which 50% of the stellar mass is formed by z — 5 
and 80% by z = 3. Our fifth model is a model from Maras- 
ton (2005). We will discuss this model further in Section 6, 
where we will use this model to test how sensitive out results 
are to our choice of stellar evolution models. 



5.2 Evolution in the stellar mass of BCGs 

In Fig. 3, we plot the observer-frame Ks-band magnitudes 
of the BCGs in our sample against their redshifts. BCGs 
from Stott et al. (2008) and Stott et al. (2010) are plotted 
as the blue and black squares, respectively, while BCGs in 
the SpARCS and CNOC1 clusters are plotted as red circles. 
While the red circles generally land within the area covered 
by the squares, the CNOC1 and SpARCS BCGs are less 
dispersed with respect to the models than the BCGs in Stott 
et al. (2008). 

Furthermore, the BCGs in the SpARCS clusters ap- 
pear to be slightly brighter than the other BCGs in the 
high-redshift subsample, although this difference seems to 
be largely driven by a few BCGs at z ~ 0.9. The clusters 
hosting the SpARCS BCGs and the other clusters in the high 
redshifts subsample have similar masses, so the correlation 
between cluster mass and the BCG stellar mass (see Sec. 5.3) 
is not the cause for the difference. Clusters in SpARCS were 
selected as galaxy overdensities, whereas most of the other 
clusters were selected through their X-ray emission. It is 
tempting to speculate that the difference is caused by the 
way the clusters were selected. However, the high-redshift 
sample is small, and we believe that a larger independent 
sample is required before one could conclusively state that 
sample selection is the reason for the difference. 

In addition to the individual BCGs, we also plot the 
predictions of the models described in the previous section. 
The normalisation 10 of the models is constrained by the data 
in the low-redshift subsample. With this normalisation, the 
BCGs at low redshift correspond to galaxies that are about 
2 mag brighter than a L* galaxy in the Coma cluster (de 
Propris et al. 1998). 

The stellar mass of individual galaxies is derived by con- 
verting the offset in magnitude between model 3, the model 
that best describes the evolution in the J-Ks colour with 
redshift, and the observed Ks-band magnitude to a stellar 
mass. We normalise the stellar mass of the BCGs to the 
stellar mass they would have by today, using the modelled 
decrease in stellar mass with time from stellar winds and 



The normalisation is computed by matching the magnitude of 
the models at the median redshift of the low-redshift subsample 
with the median magnitude of the low-redshift subsample. It dif- 
fers from the normalisation adopted in Collins et al. (2009) and 
Stott et al. (2010). In these studies the normalisation occurs over 
a more restrictive redshift interval (z < 0.05). See Sec. 6.3 for 
further details. 
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Figure 2. The observer-frame J-Ks colour of BCGs in our sample as a function of redshift. The data from this paper are plotted as 
red circles. The two most distant clusters and clusters from the CNOC1 sample are not plotted as they lack J-band data, The vertical 
dashed lines mark the boundaries of the low, intermediate and high-rcdshift subsamples that are described in the text. The evolution in 
the J-Ks colour for several stellar population models are plotted as the continuous lines. A broad range of models is shown. Note how 
well the our best fit model, model 3, which is the model we use to estimate stellar masses, describes the change in J-Ks colour with 
redshift. Sec text and Table 6 for additional details. 



Table 6. Model parameters 



Model 


Description 


Origin 


IMF 


Formation redshift 
[Gyr] 


T 


Metallicity 


1 


Low-redshift burst 


BC03 a 


Chabricr 


2.0 


0.0 


0.02 (solar) 


2 


High-rcdshift burst 


BC03 


Chabricr 


5.0 


0.0 


0.05 


3 


Best fit model 


BC03 


Chabrier 


5.0 


0.9 


60/40 split between 0.02 and 0.05 


i 


DeLucia et al (2007) model 


BC03 


Chabrier 


10.6 


1.34 


0.02 


5 


High-rcdshift burst 


M05 6 


Salpeter 


1.0 


0.0 


0.04 



Note a: Bruzual & Chariot (2003) 
Note b: Maraston (2005) 



supernova explosions (Bruzual & Chariot 2003). The results 
for the model that best follows the general evolution in the 
J-Ks colour, model 3, are shown in Table 7. In the last col- 
umn of this table, we also list the median stellar mass of 
the BCGs at the cluster redshift. When comparing the stel- 
lar mass of BCGs at low and high redshift we compare the 
stellar masses they would have by today - the second last 
column in this table. 

Without making any correction for the positive corre- 
lation between the stellar mass of the BCG and the mass 



of cluster (see Edge 1991; Burke, Collins, & Mann 2000; 
Brough et al. 2008; Whiley et al. 2008, and the next section), 
the data indicate that the stellar mass of the BCGs increase 
by a factor of 1.55 ± 0.18 between z ~ 1 and z ~ 0.17. 
The errors are determined by bootstrap resampling. The in- 
crease is found for the high-redshift subsample as a whole 
and for a smaller subsample consisting of just the clusters 
from SpARCS. In the next section we examine how the cor- 
relation between the stellar mass of the BCG and the mass 
of cluster affect our results. 
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Figure 3. The observer-frame Ks-band magnitude of BCGs as a function of redshift. The data from this paper are plotted as the red 
circles. Red circles beyond z ~ 0.8 are BCGs in the SpARCS clusters, while those below z ~ 0.8 are BCGs in the CNOC1 clusters. 
The vertical dashed lines mark the boundaries of the low, intermediate and high-redshift subsamplcs that are described in the text. 
The predicted Ks magnitudes of the models plotted in Fig. 2 and listed in Table 6 are shown as the continuous lines. The models are 
normalised to the data in the low-redshift bin. They are discussed in Sec. 5.2. Note how all models tend to underpredict the flux in 
high-redshift BCGs. 



Table 7. The three subsamples described in Sec. 4.2 and the SpARCS sample. For all quantities, we report the median value. 



Subsample 


Redshift range 


Size 


Redshift 


Cluster mass" 


Cluster Mass 6 


BCG mass c 


BCG mass d 










[10 15 M ] 


[10 15 M Q ] 


[10 12 M Q ] 


[1O 12 M ] 


Low 


z < 0.3 


93 (90) 


0.17 


0.79 e 


0.59 e 


0.45 


0.46 


Intermediate 


0.3 < z ^ 0.8 


25 (18) 


0.45 


2.34 e 


1.27 e 


0.50 


0.52 


High 


0.8 < z ^ 1.7 


32 (32) 


1.00 


1.20 


0.30 


0.29 


0.32 


SpARCS 


0.8 < 2 < 1.7 


12 (12) 


1.10 


1.19 


0.29 


0.31 


0.34 



Note a: Cluster masses corrected for the growth they are likely to have by today. 
Note b: Cluster mass at the redshift of the cluster. 

Note c: The stellar mass of the BCG at redshift zero (accounts for stellar mass loss). 
Note d: The stellar mass of the BCG at the redshift of the cluster. 

Note e: Computed for the subset of clusters (numbered in brackets) with masses. See text for details. 



5.3 Accounting for cluster masses 



The stellar mass of BCGs correlates with cluster mass in the 
sense that larger clusters tend to host larger BCGs (Edge 
1991; Burke, Collins, & Mann 2000; Brough et al. 2008; 
Whiley et al. 2008; Stott et al. 2012). Comparing the stellar 
mass of BCGs in our three subsamples without accounting 



for this correlation will lead to biased results if the median 
mass of the clusters in the subsamples differ. 

To account for this correlation, we first need to estimate 
how clusters grow in mass so that we can fairly compare clus- 
ters that are observed at different redshifts. Over the redshift 
range that our sample covers, clusters grow significantly. For 
example, in the hierarchical model of structure formation, 
a cluster with a mass of 5 x 10 14 Mq at z = 1 is predicted 
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to grow by a factor of about three by today (Wechsler et 
al. 2002; Fakhouri, Ma, & Boylan-Kolchin 2010). Fakhouri, 
Ma, & Boylan-Kolchin (2010), who use the Millennium and 
Millennium II simulations (Springel et al. 2005; Boylan- 
Kolchin et al. 2009) and Wechsler et al. (2002), who use an 
independent simulation (Bullock et al. 2001), find similar 
growth rates. We use mean accretion rates in Fakhouri, Ma, 
& Boylan-Kolchin (2010) to estimate the mass each cluster 
should have by the current epoch using the masses they had 
at the redshifts they were observed. We describe how we 
estimate cluster masses at the redshift at which they were 
observed in Sec. 5.3.1. 

After evolving our clusters forward in time to today, we 
find that the median mass of the clusters in our three sub- 
samples differ by as much as a factor of three (see Table 7). 
The differences in the subsamples reflect the volumes probed 
and sensitivity limits of the surveys that were used to build 
our subsamples. Since clusters in the intermediate and high- 
redshift subsamples are, by the current epoch, more massive 
than those in the low-redshift subsample, the correlation be- 
tween cluster mass and BCG stellar mass - if uncorrected - 
leads to an underestimate in the amount of evolution in the 
stellar mass of BCGs. 

In this paper, we explore a couple of approaches to 
account for the correlation. In the first approach, we first 
match the cluster mass distributions in the samples being 
compared before comparing the masses of the BCGs. In the 
second approach, we normalize BCG stellar masses to some 
fiducial mass using the relationship between cluster mass 
and the BCG stellar mass. 

Our approaches to account for this correlation differs 
from approaches used in the past. In Whiley et al. (2008), 
clusters are grouped according to the mass they had at the 
redshift they were observed. In Stott et al. (2010), the mass 
of BCGs are compared to the mass of BCGs from semi- 
analytic models after first matching the masses of the clus- 
ters in the semi-analytic models to the observed masses. 



5.3.1 Estimating Cluster mass 

Given the heterogeneous nature of the data that is available 
for our clusters, we estimate cluster masses 11 in different 
ways. For clusters in the low and intermediate-redshift sub- 
samples of Stott et al. (2008) , we use the M500 masses listed 
in Mantz et al. (2010). These masses are converted to M200 
assuming that the cluster mass profile follows a Navarro- 
Frenk- White profile (Navarro, Frenk, & White 1997) with a 
concentration index of 5. For the conversion, we use the for- 
mulae listed in Appendix C of Hu & Kravtsov (2003). Only 
about a quarter of the clusters in Stott et al. (2008) are in 
Mantz et al. (2010). To increase the number of clusters in 
Stott et al. (2008) with mass measurements, we use the X- 
ray temperatures (Ebeling et al. 2007, 2010) and the X-ray 
luminosities (Ebeling et al. 1996, 1998, 2000) of these clus- 
ters (the luminosities were corrected for the cosmology used 



11 We use M200 for cluster masses. M200 is the mass contained 
within a radius within which the mean density of the cluster ex- 
ceeds the critical density of the Universe at the redshift of the 
cluster by a factor of 200. 



in this paper) and the temperature-mass and luminosity- 
mass relations in (Mantz et al. 2010) to estimate cluster 
masses. 

Twenty one of the clusters in (Ebeling et al. 2007, 2010; 
Mantz et al. 2010) are also in (Ebeling et al. 1996, 1998, 
2000), which enables us to compare the mass derived from 
the X-ray luminosity with the mass derived from X-ray tem- 
perature. The median ratio is 1.14 with a s.d. of 0.53. We 
correct the masses determined from the X-ray luminosity 
by the median ratio and use the s.d. as a measure of the 
uncertainty in the conversion. 

For the SpARCS clusters, we estimate the cluster mass 
from the line-of-sight velocity dispersion (Wilson et al., in 
preparation). For clusters in the high-redshift subsample of 
Stott et al. (2010), we use the X-ray temperature reported 
in that paper and convert them to masses using the relation 
in Mantz et al. (2010). For all but two of the clusters in 
the intermediate-redshift subsample of Muzzin et al. (2007a) 
we use the X-ray temperatures listed in Hicks et al. (2006) 
and convert them to masses using the relation in Mantz 
ct al. (2010). For the remaining two clusters (MS1224+20 
and MS1231+15), we use the masses listed in Muzzin et al. 
(2007b) , which are computed from the line-of-sight velocity 
dispersion. 

After removing clusters from Stott et al. (2008) that are 
not listed in Ebeling et al. (1996, 1998, 2000, 2007, 2010) and 
Mantz et al. (2010), we end up with 90 and 18 clusters in 
the low- and intermediate-redshift subsamples, respectively. 
The number of clusters in the high-redshift subsample is 
unchanged. The numbers are listed in Table 7. Altogether, 
there are 152 clusters in our three subsamples. 

5.3.2 Cluster mass vs. BCG stellar mass 

The correlation between cluster mass (at the redshift at 
which it was observed) and BCG stellar mass for these three 
subsamples is shown in Fig. 4. Errors in the mass of the 
BCGs are derived from errors in the photometry. If an er- 
ror in the Ks-band photometry was unavailable, we con- 
servatively set the error to 10%. Our results are not very 
sensitive to this value, as errors in the cluster masses are 
much larger. Errors in the mass of the clusters are discussed 
in 5.3.3. We note that clusters in the intermediate-redshift 
subsample generally have higher masses than clusters in the 
low-redshift subsample. As noted earlier, the difference be- 
tween subsamples reflect the volumes probed and the sen- 
sitivity limits of the individual surveys that were used to 
build the subsamples. 

We fit a power law to the data, using a lognormal distri- 
bution to represent the likelihood of getting a certain data 
point given the model and allowing for additional dispersion 
by scaling the measurement uncertainties. The index of the 
power law that corresponds to the maximum of the poste- 
rior distribution is 1.6 ±0.2. Because we treat errors and the 
amount of extra dispersion in both axes equally, our results 
are robust to flipping the axes in the fit. The amount of ex- 
tra dispersion found in the fit corresponds to increasing the 
size of the error bars by a factor of 1.5. 

The index of the power law suggests that clusters ac- 
crete mass five times faster the BCGs accrete stellar mass. 
Within uncertainties, the index is similar to that found in 
Stott et al. (2010), who find 2.4±0.6 and Stott et al. (2012) 
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Figure 4. The correlation between the mass of the cluster at 
the epoch at which it was observed and the stellar mass of the 
BCG. The different symbols represent different redshift ranges. 
The solid line is a fit to the data. Setting the index of the power 
law to the value reported in Hansen et al. (2009) results in a 
poorer fit to the data (dashed line). 



who find 1.3±0.1. Some of the difference between our results 
and those in Stott et al. (2010, 2012) come from the way the 
samples are selected and the way the analysis is performed. 
Our best fit index is about a factor of two smaller than those 
reported in earlier works (Lin & Mohr 2004; Popesso et al. 
2007; Brough et al. 2008; Hansen et al. 2009). For exam- 
ple, Hansen et al. (2009), find an index of 3.3 between the 
i-band luminosity (k-corrected to z = 0.25) and M200. We 
redid the fit with the index constrained to this value. The 
resulting relation is shown in Fig. 4 as the dashed line. It is 
a poorer fit to the data. 

Our fit to the entire sample seems to be largely driven 
by the clusters in the high-redshift subsample, whereas most 
of the clusters in Hansen et al. (2009) were at low redshift. 
This raises the possibility that there is evolution in the index 
of the power law with redshift. Alternatively, the difference 
might be caused by redshift-dependent selection effects. In 
accounting for cluster masses in the following sections, we 
adopt a conservative approach and examine how our results 
depend on which index we choose to use. We will find that 
our conclusions are robust to this choice. 



5.3.3 Accounting for cluster masses 

As foreshadowed earlier, we use two approaches to account 
for the correlation between cluster mass and BCG stellar 
mass. We discuss the first approach in this section and dis- 
cuss the second approach in the section that follows. 

In the first approach, we randomly select clusters from 
the three subsamples until the mass histograms 12 of the sub- 
samples match. Clusters are matched according to the mass 
they will have by the current epoch. Implicit in this ap- 
proach is the method we use to estimate how clusters build 
up their mass with time. 

We cannot match all three subsamples simultaneously, 

12 We use a bin size of 2 X 10 14 and we use the mass the clusters 
are likely to have by today. 



Figure 5. Upper Panel A histogram showing the distributions 
of cluster mass (extrapolated to the current epoch) for the low 
and high-redshift subsamples. The median masses are marked 
with the downward pointing arrows. Note how the median mass 
of the two distributions differ and how skewed the low redshift 
subsample is with respect to the high-redshift one. Lower two 
panels Histograms of the re-sampled low and high-redshift sub- 
samples. They are resampled so that they are identical for a bin 
width of 2 X 10 14 Mq . The median masses, marked with the down- 
ward pointing arrows, are now more similar. There are 23 objects 
in the lower two histograms. 



because trying to get all the histograms to match would 
result in very few objects per subsample. Instead we compare 
the low-redshift subsample with the intermediate and high- 
redshift subsamples separately. The method is illustrated in 
Fig. 5 for the comparison between the low and high-redshift 
subsamples. 

In order to get a measure of the uncertainties in the 
derived mass ratios, we do two things. We first perturb the 
cluster mass by an amount that depends on two sources of 
error: the uncertainty in the measurement of the mass proxy 
(X-ray temperature, X-ray gas mass, X-ray luminosity or 
line-of-sight velocity dispersion) and the intrinsic scatter 
between the mass proxy and the mass. For masses that are 
determined from the X-ray gas mass or the X-ray temper- 
ature, we assign a scatter of 15% (Mantz et al. 2010). For 
masses determined from the line-of-sight velocity dispersion, 
we assign a scatter of 30% (Hicks et al. 2006). For masses 
inferred from the X-ray luminosity, we use 50%, which we 
derived earlier. The magnitude of the perturbation is drawn 
from a lognormal distribution. The s.d. of the distribution 
is set equal to the two uncertainties added in quadrature. 

Secondly, we resample the three subsamples with re- 
placement (bootstrap resampling) to allow for uncertainties 
that come from sample size. Only then do we try to match 
the histograms in the three subsamples. We repeat this exer- 
cise 100 times for each comparison to create 200 realisations 
from the data. For each realisation, we compute the median 
BCG stellar mass, the median cluster mass and median red- 
shift. For each comparison, we then average the results from 
the 100 realisations and get an estimate of the robustness of 
the results from the variance. The results of the comparisons 
are listed in Table 8 and shown in Fig. 6. The uncertainty in 
the last column in Table 8 is computed from the 100 realisa- 
tions and gives an indication of the robustness of the result. 



© 2002 RAS, MNRAS 000, 1-?? 



12 C. Lidman et al. 



▲ Matched histograms 




'8.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 
redshift 



Figure 6. The evolution in the median stellar mass of BCGs as 
a function of redshift. The green triangles take into account the 
correlation between cluster mass and the stellar mass of its BCG 
by matching clusters according to the masses they will have by 
the present epoch. In a second approach, the red and blue circles 
account for this correlation using the relations shown in Fig. 4. 
The small black squares do not account for this correlation. Note 
how all the points in the high-redshift bin lie below the red line, 
how the green, blue and red points in the intermediate and high- 
redshift bins lie below the black points, and how these points are 
a better match to the De Lucia &; Blaizot (2007) model (solid 
blue line). All points are normalised so that their low-redshift 
points land on this model. The vertical dashed lines mark the 
boundaries of the low, intermediate and high-redshift subsamples 
that are described in the text. The points are plotted at the me- 
dian redshifts of the subsamples. They differ slightly between the 
green, red, blue, and black points because a more restricted range 
of clusters is selected when matching cluster masses. See text for 
details on how the error bars are computed. The red horizontal 
line represents no mass evolution. The data used in this plot are 
summarised in Table 10 

The uncertainties are plotted as the vertical error bars in 
Fig. 6. No other uncertainties are included in these error 
bars. 

Between z ~ 0.9 and z ~ 0.2, the stellar mass of BCGs 
increase by a factor of 1.8 ± 0.3. This is larger than the 
increase reported in the previous section, which did not ac- 
count for the correlation between cluster mass and the stellar 
mass of the BCG. 

We repeated our analysis by comparing the stellar mass 
of BCGs in clusters that have the same mass at the redshift 
they were observed. This is the method used in Whiley et 
al. (2008). The results are presented in Table 9. Not sur- 
prisingly, due to the correlation between the stellar mass of 
the BCG and the mass of the cluster, and the considerable 
growth in cluster mass between z ~ 1 and today, the evo- 
lution in the stellar mass of the BCG is less evident when 
clusters are compared in this way. 

5.3.4 An alternative approach 

An alternative approach to account for the correlation be- 
tween cluster mass and BCG stellar mass is to adjust the 
BCG stellar mass according to the relation shown as the 
solid line in Fig. 4. As in the first approach, we use the clus- 
ter mass extrapolated to the current epoch and not the mass 



Table 10. A summary of the data appearing in Fig. 6 



Method 


Redshift 


Mass Ratio" 


Matched histograms 


0.20 


0.81 




0.40 


0.77 ± 0.12 




0.96 


0.47 ± 0.07 


No matching 


0.17 


0.84 




0.45 


0.97 ± 0.12 




1.00 


0.54 ± 0.06 


Correcting with index 1.67 


0.17 


0.84 




0.45 


0.48 ± 0.10 




1.00 


0.36 ± 0.06 


Correcting with index 3.33 


0.17 


0.84 




0.45 


0.73 ± 0.10 




1.00 


0.47 ± 0.06 



Note a: The ratios are scaled so that the low-redshift point 
matches the prediction of De Lucia &; Blaizot (2007). 



they had at the epoch they were observed. The results are 
shown as the red points in Fig. 6. The errors are derived 
using bootstrap resampling. 

Compared to the previous approach, we find stronger 
growth in the BCG stellar mass as a function of redshift and 
better agreement between the data and the semi-analytic 
models of De Lucia & Blaizot (2007). However, the result 
depends on the relation shown in Fig. 4. Adopting the re- 
lation found by Hansen et al. (2009) instead of the relation 
we find, for example, results in less growth (the blue points 
in Fig. 6). The point at intermediate redshifts is affected 
most, since these clusters will be, by today, three times more 
massive than clusters in the low-redshift sample (see Ta- 
ble 7), thereby leading to significant adjustments. The point 
at high-redshifts is affected less, because these clusters will 
be, by the current epoch, similar in mass to clusters in the 
low-redshift sample. 

We do not adopt the relation found in Hansen et al. 
(2009). Instead we use it to demonstrate the sensitivity of 
the approach to changes in the power-law index. The sam- 
ple used in Hansen et al. (2009) to compute the relation is 
restricted to clusters in the redshift range 0.1 ^ z ^ 0.3. 
There are also differences in the analysis. Cluster masses in 
Hansen et al. (2009) are estimated from the optical richness 
and BCG masses are estimated from the observer-frame i— 
band, which, when combined with the optical selection, may 
lead to biases that influence the result. 

We use the differences in the results between this ap- 
proach (using the best fit power law index), the approach 
described in Sec. 5.3.3, and the approach of not applying any 
correction as an estimate of the size of the systematic error. 
Clearly, between the intermediate and low-redshift samples, 
the evidence for evolution is marginal. The statistical and 
systematic uncertainties are too large. 

However, between the high and low-redshift samples, 
the evidence for evolution is clear and unambiguous. Be- 
tween z ~ 0.9 and 2 ~ 0.2, the stellar mass of BCGs increase 
by a factor of 1.8 ± 0.3 with a spread of 0.4 spanning the 
three approaches. 



© 2002 RAS, MNRAS 000, 1-?? 



Brightest Cluster Galaxies 13 

Table 8. BCG mass ratios. The cluster mass distributions are matched using the masses the clusters will have by the current epoch. 



Samples a and b 


Median Rcdshifts 


Median cluster masses 


Median BCG masses 


BCG mass ratio 1 




Sample a 


Sample b 


Sample a Sample b 


Sample a Sample b 










[10 15 M ] 


[1O 12 M ] 




low (a) - intermediate (b) 


0.20 


0.40 


1.59 1.61 


0.54 0.51 


0.96 ± 0.20 


low (a) - high (b) 


0.17 


0.97 


0.84 0.84 


0.45 0.26 


0.58 ± 0.08 



Note 1: Defined as the median stellar mass of subsample b divided by the median stellar mass of subsample a 



Table 9. As for Table 8 with the difference that the matching is done using the masses the clusters have at the redshift at which they 
were observed. 



Samples 


a and b 


Median Rcdshifts 


Median cluster masses 


Median BCG masses 


BCG mass ratio 1 






Sample a 


Sample b 


Sample a Sample b 


Sample a Sample b 












[10 15 M Q ] 


[10 12 M ] 




low (a) - 


intermediate (b) 


0.20 


0.44 


1.02 1.06 


0.51 0.51 


1.02 ± 0.16 


low (a) - 


- high (b) 


0.10 


1.00 


0.33 0.30 


0.41 0.29 


0.71 ± 0.10 



Note 1: Defined as the median stellar mass of subsample b divided by the median stellar mass of subsample a 



6 DISCUSSION 

The semi- analytic model of De Lucia & Blaizot (2007) pre- 
dicts that BCGs grow by a factor of almost three in stellar 
mass between z ~ 0.9 and z ~ 0.2. Over the same redshift 
interval, we observe that BCGs increase their stellar mass 
by a factor of 1.8 ± 0.3. 

Our result depends on the methods we have used to 
analyse the data and the choice of models that we have used 
to estimate stellar masses. We discuss each of these in turn, 
finding that our results are robust. 



6.1 Estimating the Ks band flux 

Throughout this paper we have estimated the Ks-band 
flux using MAG_AUT0 in SExtractor, which is a Kron-likc 
magnitude (Kron 1980) within an elliptical aperture. Un- 
doubtably, MAG_AUT0, like all other measures of the total 
magnitude used in the literature, will be systematically bi- 
ased to low or high values depending on the nature of the 
object being measured. The bias can occur for a number of 
reasons, such as the number and brightness nearby neigh- 
bours, the presence of intra-cluster light and/or a cD en- 
velope, the point spread function (i.e. seeing) and residual 
errors that come from the imprecise removal of the bright 
night sky from near-IR images. For the purpose of compar- 
ing the stellar mass of BCGs at low and high redshifts, the 
most important aspect of the measurement is that the bias 
does not change with redshift. In this section we search for 
evidence of the bias changing with redshift and try to esti- 
mate how large this bias may be. 

Graham & Driver (2005) noted that Kron-like magni- 
tudes can significantly underestimate the flux of galaxies 
with Sersic profiles. We ran our own investigation into the 
accuracy of MAG_AUT0 by inserting objects with Sersic pro- 
files in simulated images. In the simulations, we mimicked 
the background noise and image quality of the real data. 



For de Vaucouleur profiles 13 , MAG_AUT0 misses between 18% 
and 35% of the flux, depending on the redshift. The trend 
with redshift is non- monotonia At z = 0.1, 25% of the flux 
is missed. This decreases to 18% by z = 0.4, then increases 
to 25% by z = 1.0 and to 35% by z = 1.6. For higher Sersic 
indices (we tested indices as high as n = 8), higher fractions 
of the flux are missed by MAG_AUT0; however the trend with 
redshift is the same. 

Clearly, if the profiles of BCGs evolve with time, then 
there will a redshift-dependent bias in the stellar masses 
that are inferred from the photometry. For example, if low- 
redshift BCGs had de Vaucouleur profiles and high-redshift 
BCGs had Sersic profiles with n = 8, then we would over- 
estimate the flux of the low redshift BCGs relative to their 
distant cousins, and therefore their stellar mass, by around 
10%. If the opposite was true (i.e. high-redshift BCGs had 
de Vaucouleur profiles and low-redshift BCGs had Sersic 
profiles with n = 8), then we would underestimate the flux 
of the low redshift BCGs by 3%. The asymmetry is caused 
by the dependence of how accurately MAG_AUT0 measures to- 
tal magnitude with redshift and seeing. Observational con- 
straints on the redshift dependence of the Sersic index show 
that the redshift dependence is much weaker than the range 
of values that we have considered here (Stott et al. 2011). 

In our simulations, we neglected errors in the photome- 
try that come from nearby (in projection) galaxies and im- 
precise sky-subtraction. To investigate these issues, we com- 
pare the integrated light profiles of the BCGs measured with 
version 3.0.4 of GALFIT 14 (Peng 2002). Galaxies neighbour- 
ing the BCG were either fitted simultaneously (if they were 
within 2-3" of the BCG) or masked as bad pixels (if they 
were further than this). We used the residual images and the 
reduced \ 2 °f the fit to determine how well the data were 
described with the model. With the exception of SpARCS- 
0035 and SpARCS-1638, the BCGs could be modelled sat- 
isfactorily (a reduced \ 2 close to one) with a de Vaucouleur 



13 The de Vaucouleur profile is equivalent to a Sersic profile with 
the Sersic index, n, set to 4. 

14 http://users.obs.carnegiescience.edu/peng/work/galfit/galfit.html 
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profile. The BCGs of both SpARCS-0035 and SpARCS-1638 
were better fit with Sersic profiles that had a higher Sersic 
indices. 

For the SpARCS clusters, which have a median redshift 
of 2 ~ 1.1, the offset between MAG_AUT0 and the magnitude 
determined by integrating the fitted GALF1T profile out to 
infinity has a median value of 0.30 mag (i.e. relative to the 
integrated GALFIT flux, MAG_AUT0 underestimates the flux. 
For clusters in the CNOC1 sample, which have a median red- 
shift of 2 ~ 0.28, the offset between the integrated GALFIT 
magnitude and MAG_AUT0 has a median value of 0.49 mag, 
which is considerably larger than the median value found 
for the SpARCS clusters. 

In part, the difference between the offsets is due to the 
way the aperture in MAG_AUT0 is defined. The size of the aper- 
ture depends on the seeing convolved profile of the BCG. 
Since the CNOC1 and SpARCS samples were taken in simi- 
lar seeing, the apertures for the SpARCS BCGs are affected 
more by the seeing, since the angular size of the BCGs rel- 
ative to the seeing disk is smaller. This leads one to using 
apertures for the SpARCS BCGs that are larger than one 
would have used if the ratio of the seeing to the angular size 
of the BCG was the same for both samples. This then leads 
to a smaller difference for SpARCS BCGs. 

The offset between the differences translates directly 
into a relative offset in the stellar masses of the BCGs in 
SpARCS and CNOC1. Relative to the stellar masses of the 
BCGs in the SpARCS clusters, we are underestimating the 
stellar masses of the BCGs in the CNOC1 clusters by a 
factor of 1.2. If this offset were applicable to the rest of 
the BCGs in our low-redshift subsample, then we would be 
underestimating the growth in BCGs between the low and 
high-redshift subsamples by a similar amount. Hence, in- 
stead of finding that the mass grows by a factor of 1.8 ± 0.3, 
we would find that the mass grows by a factor of 2.2. 

In this paper, we do not use the integrated magnitude 
in GALFIT to estimate stellar masses. We make this choice 
because most of the BCGs in our low-redshift subsample 
have not been analysed with GALFIT. 

Instead, we note that there is a source of systematic 
uncertainty in the relative stellar masses between low and 
high-redshifts that comes from the photometry. By compar- 
ing two widely used techniques to do galaxy photometry, we 
estimate this uncertainty to be ~ 20%. 

6.2 Stellar masses of the BCGs 

In section 5.2, we described how we used the Ks-band mag- 
nitude of BCGs and the predictions of a model that broadly 
describes the change in the J-Ks colour of BCGs with red- 
shift to estimate their stellar masses. The masses will depend 
on the model used, so our conclusions are model dependent. 
To explore how sensitive this dependence is, we re-estimate 
the masses using another stellar population synthesis code. 

Model 5 in Fig. 2 and Table 6 is from Maraston (2005, 
hereafter M05). The stars in this model formed in a single 
burst at z = 4, have a metallicity that is twice solar and form 
a red horizontal branch, as found in most metal-rich globular 
clusters. The model follows the evolution of the J—Ks colour 
with redshift almost as well as model 3, although it tends to 
predict redder colours at z ~ 1.4. 

The M05 and BC03 models differ in several ways. One 



of the differences most noted in the literature (see Maraston 
et al. 2006; Marigo et al. 2008, for example) is the treatment 
of the thermally pulsing asymptotic giant branch (TP-AGB) 
phase of stellar evolution. This phase significantly affects 
the optical-NIR colours of simple stellar populations in the 
age range 0.5 < t < 1.5 Gyr. Over this age range, the M05 
models predict redder optical-NIR colours and, for a given 
stellar mass, higher NIR luminosities. 

As we did for the BC03 models, we normalise the M05 
model so that they match the brightness of the BCGs over 
the redshift interval < z < 0.3. Using the M05 models, we 
find that the stellar mass of BCGs at z = 0.9 are 1.81 ± 0.26 
times less massive than BCGs at z = 0.2, which is similar 
to the results that we derive using the BC03 models. 

We do not know if the BC03 models are more appro- 
priate than those in M05; however, we note that recent ob- 
servations are now suggesting that the contribution from 
TP-AGB stars to the near-IR flux, may not be as significant 
as previously thought (Kriek et al. 2010; Zibetti et al. 2012). 

Finally, for a couple of BCGs, we examine how well 
our stellar mass estimates compare with measurements that 
are made using more extensive photometric data, a differ- 
ent stellar population model and a different way of estimat- 
ing total magnitudes. Using 10 broad band filters extending 
from 4640 A to 8.0 {im (rest frame), Rettura et al. (in prepa- 
ration) estimate a stellar mass of 3.9 x lO n M for the BCG 
in SpARCS- J003550-431224. Our estimate from the K band 
photometry is 3.3 x lO n M . Rettura et al. (2006), using 
9-band photometry, derive a stellar mass of ~ 2.3 x lO n M 
for the BCG in RDCS J1252. 9-2927. From the K band pho- 
tometry, we find 3.3 x lO n M . 

6.3 Comparison with other results 

Our finding of significant evolution in the stellar mass of 
BCGs with time differs from the findings of a number of 
authors (Whiley et al. 2008; Collins et al. 2009; Stott et al. 
2010). Since the Ks-band magnitudes and J-Ks colours of the 
BCGs that we have added in this paper are similar to the 
magnitudes and colours of BCGs from earlier works (Stott 
et al. 2008, 2010) and since much of our sample consists of 
BCGs from these works, the reason for the difference lies in 
the way we have done the analysis. 

In part, the difference comes from the way we have 
compared low and high-redshift BCGs. In this paper, we 
first match clusters according to the mass they will have 
by the current epoch, before comparing the stellar mass of 
the BCGs they host. Earlier works have done this compari- 
son differently. For example, Whiley et al. (2008) match the 
clusters according to the mass they had at the redshift they 
were observed. We repeated our analysis using this approach 
(see Table 9 and Sec. 5.3.3) and found that the evidence for 
evolution became considerably weaker. 

The difference may also come from the redshift inter- 
val that we use to define the low-redshift subsample. In this 
paper, we use z < 0.3. This is broader than that used by 
other authors, e.g. Stott et al. (2008). The broader inter- 
val allows us to use more objects to determine the stellar 
mass ratio at the expense of a smaller time interval between 
the low and high-redshift samples. We repeated our analysis 
with the redshift interval for the low-redshift sample set to 
2 < 0.1. We find that the stellar mass ratio between the high 
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and low-rcdshift subsamples increases slightly to 0.61 ±0.19. 
The uncertainty is larger because there are fewer objects 
in the low-redshift subsample. More significant, however, is 
that the redshift interval between the low and high-redshift 
subsamples increase, thereby increasing the tension between 
the data and the predictions of the semi-anlaytic models. 

6.4 The build up of stellar mass in BCGs 

In semi-analytic models of (De Lucia & Blaizot 2007), the 
stellar masses of BCGs increase by a factor of about 3 be- 
tween z ~ 0.9 and z ~ 0.2. Our results suggest that the 
growth is slower than this. Over the same redshift range, 
we find the increase to be a factor of 1.8 ± 0.3, suggesting 
that the model over-predicts the amount of stellar mass by 
a factor of ~ 1.5 - the difference between the blue line and 
the green triangle in Fig. 6). 

In semi-analytic models, most of the build up in stellar 
mass occurs through dry mergers (both major and minor) 
with other galaxies. There is ample observational evidence 
for major mergers in the centres of clusters (Rasmussen et 
al. 2010; Brough et al. 2011; Bildfell et al. 2012). Brough et 
al. (2011) in a study of three BCGs at z ~ 0.1 with nearby 
companions found that the companions of two of the BCGs 
would merge with the BCG within 0.35 Gyr. More dramatic 
still is the the merger that is occurring in the centre of MZ 
10451 (Rasmussen et al. 2010). 

Evidence for major mergers can also be seen in some 
of the SpARCS BCGs, For example, the isophotes of the 
BCG in SpARCS J163435+402151 are distorted, indicating 
a possible major merger. In this cluster, there is evidence 
that another major merger is occurring for a galaxy that is 
almost as bright on the other side of the cluster. 

There is also an example in the SpARCS sample of a 
merger that is likely to happen by today. In the centre of 
SpARCS J161641+554513, there is a galaxy that is within 
20kpc projected distance of the BCG. The velocity differ- 
ence between the two galaxies is ~ 140km/s, and the com- 
panion is almost as bright as the BCG (see Fig. 1). It is 
highly likely that these two galaxies would have merged 
by now. Additional examples of likely major mergers can 
be found in RX J0848.9+4452 (Yamada et al. 2002) and 
RDCS J1252-2927 (Collins et al. 2009). 

While it is clear that mergers do occur, it is not yet clear 
what fraction of the stars in the merging galaxies end up in 
the BCG and what fraction end up distributed throughout 
the cluster, appearing as intra-cluster light (ICL). The ap- 
parent lack of evolution in the stellar mass of BCGs that 
was found in earlier work suggested that the contribution 
to the ICL was close to 100% (Whiley et al. 2008; Collins 
et al. 2009; Stott et al. 2010). Our results suggest that it is 
closer to 50%. High resolution simulations suggest that 50- 
80% of the mass of mergers will be distributed throughout 
the cluster (Conroy, Wechsler, & Kravtsov 2007; Puchwein 
et al. 2010). Recent measurements of the ICL show that the 
ICL grows relative to the total cluster light by a factor of 
2-4 since z ~ 1 (Burke et al. 2012). 

It is possible that some of the BCGs in our sample 
are increasing their stellar mass through star formation. 
Out of the 12 SpARCS BCGs, 5 show emission from the 
[Oil] AA 3726,3728 doublet, which is an indicator of star for- 
mation and/or AGN activity. It is unlikely that most of the 



[Oil] emission that we detect comes from star formation. 
Over 70% of low redshift BCGs that have detectable [Oil] 
emission have line ratios that are consistent with the line ra- 
tios of AGN (von der Linden et al. 2007). Our spectra do not 
cover the lines that can be used to separate between AGN 
activity and star formation, such as the [OIII] AA 4959,5007 
doublet, H/3, Ha and [NIIJA6584. 

If we were to assume that the [Oil] emission did come 
from star formation entirely and if we ignore dust, then the 
average [Oil] line flux corresponds to a star formation rate 
of about 1 solar mass per year, using the equation (4) in 
Kewley, Geller, & Jansen (2004) to make the conversion be- 
tween [Oil] line flux and the star formation rate. At these 
rates, star-formation will not contribute much to the over- 
all stellar mass of the BCG, even if they were to continue 
forming stars at this rate until today. 



7 SUMMARY AND CONCLUSIONS 

Using near-IR photometry from the literature (Stott et al. 
2008, 2010) and photometry from an analysis of imaging 
data that we obtained using several ground-based near-IR 
cameras, we have investigated how the stellar masses of 
BCGs change with redshift. The BCGs in our sample cover 
a broad redshift range, from z — 0.03 to z — 1.63, which 
covers 9.8 Gyr, or 70% of the history of the universe. 

To estimate the stellar mass of the BCGs, we compare 
the Ks band flux with the predictions from a stellar popu- 
lation synthesis model that matches the J-Ks colour of the 
BCGs over the entire redshift range covered by the data. 

We then compare mass of BCGs at low and high red- 
shifts. After accounting for the correlation between BCG 
stellar mass and cluster mass, we find that, between z = 0.9 
and z = 0.2, BCGs, on average, grow in mass by a fac- 
tor of 1.8 ± 0.3. Our result is not weakened if we choose 
other methods to estimate the Ks band flux or if we choose 
other stellar population synthesis models to infer the mass. 
The systematic uncertainty coming from the photometry is 
probably the dominant source of systematic uncertainty in 
our analysis and affects our estimates of the growth rate by 
around 20%. 

Our conclusions differ from those of earlier works 
(Collins et al. 2009; Stott ct al. 2010). In part, this is due to 
the way we have accounted for the correlation between the 
mass of the BCG and the mass of the cluster and to the red- 
shift intervals that we use to define the low and high-redshift 
subsamples. 

Our measurements are now in better agreement with 
the predictions of semi-analytic models for the growth of 
stellar mass in BCGs (De Lucia & Blaizot 2007). However, 
there is still some tension between the data and these mod- 
els, which predict growth rates that are a factor of 1.5 higher. 

We find direct evidence that some of the BCGs in 
our sample are building up their stellar mass through 
star-formation and major mergers. However, star-formation, 
while present in some of BCGs, is at low levels. At these lev- 
els, star formation cannot be the dominant mechanism for 
the build up of stellar mass in BCGs over the last 10 bil- 
lion years. The build-up mainly occurs through mergers, of 
which some are clearly major. 
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APPENDIX A: DATA 

In Tables Al, A2 and A3, we list the names of the clusters, 
their redshifts, their masses and the mass proxy used to de- 
termine masses. Two estimates of the mass are provided. 
The first is computed from the mass proxy and represents 
the mass of the cluster when it was observed. The errors in 
the cluster masses only include the error in the mass proxy. 
They do not include the intrinsic scatter in the relation be- 
tween mass and mass proxy. The second mass extrapolates 
the first mass to the current epoch by integrating the mean 
mass accretion rates in Fakhouri, Ma, & Boylan-Kolchin 
(2010). The error does not take into account the intrinsic 
scatter in the accretion rates. Wechsler et al. (2002) esti- 
mates that between 2 = 1 and z = 0, the final mass of a 
10 14 M Q halo can scatter by 20-30%. Also listed in the table 
are the magnitudes and colours of the BCGs. If available, we 
also list the errors in these quantities. Excluding the BCGs 
in the SpARCS and CNOC1 clusters, the magnitudes and 
colours of the BCGs in these tables were obtained from Stott 
ct al. (2008, 2010) and J. P. Stott (private communication). 
The model-dependent masses, which have been adjusted to 
account for the loss of mass due to supernova explosions and 
stellar winds, are listed in the final column. Not all clusters 
have mass measurements in the papers listed in the main 
body of the paper. These clusters are listed in Table A4. 
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Table Al. The low-redshift subsamplc 







Cluster 








BCG 




Name 


Rcdshift 


Mass 


Mass today 


Mass Proxy 


Ks 


J-Ks 


Stellar Mass 






[10 15 M ] 


[10 15 M ] 






[magj 


[magj 


Tin 12 Mr-,] 


Abelll902 


0.160 


0.48i°o°06 


o.6iin 8 8 


X-ray 


luminosity 


12.63 


1.37 


0.59 


Abelll93 


0.049 


0.181°"™ 


o.2oin 3 3 


X-ray 


luminosity 


10.43 


1.06 


0.48 


Abelll930 


0.131 


0-381^5 


0.48+g-g; 


X-ray 


luminosity 


12.47 


1.10 


0.48 


Abelll991 


0.059 


16+ 002 
u - lu -0.02 


0.l8l°;° 2 


X-ray 


luminosity 


11.15 


1.01 


0.35 


Abell2029 


0.077 


i n+0.2 
— 2 


i-4l° : 2 2 


X-ray gas mass 


10.30 


1.13 


1.31 


Abell2034 


0.113 


0.89l°o« 


-1 9 +0.2 
1,z — 2 


X-ray 


ido llldoo 


12.21 


1.07 


0.46 


Abell2052 


0.035 


0.25l° ° 2 


0.28lg;° 2 


X— ray 


li i m i Tiric;i1""v^ 

lLllllllll_Ji3lLry 


9.88 


1.00 


0.41 


Abell2065 


0.073 


43+ - 04 
—0 04 


0.48±°"g| 


X— ray 


i ui iiin u&iir y 


12.03 


1.15 


0.24 


Abell2072 


0.127 


31+° : ° 6 


39+ 08 


X— ray 


li i m 1 vii~4Qi1""vr 

1L111111H.J01L y 


12.82 


1.23 


0.33 


Abell2107 


0.041 


10+0^02 
Ull ° — 0.02 


U -- LO _0 02 


X— ray 


1 LlllllllU&lLry 


10.10 


1.00 


0.45 


Abell2124 


0.066 


16+ 003 
U - 1U — 0.03 


o.i7in 3 3 


X— ray 


1 LlllllllU&lLry 


11.05 


1.05 


0.48 


Abell2175 


0.095 


29+ 03 
u-^f _q 03 


0-32±°-g* 


X— ray 


1 1 1 TY1 1 Tl(~!C;i1""V? 

i uiiiimjoiir^y 


11.78 


1.18 


0.50 


Abell2204 


0.152 


l-4l° 1 


1.8l°- 3 


Y-l'flV 


^as mass 


12.23 


1.14 


0.78 


Abell2244 


0.097 


0.821°;!* 


0.94l°-i 7 7 


X-ray 


^cxa llldSB 


11.59 


1.11 


0.62 


Abell2259 


0.164 


o ^+°- 08 


0.711°!° 


X— ray 


1 LlllllllUblLry 


12.77 


1.19 


0.54 


Abell2345 


0.177 


76+ al3 

u - ' u — 0.13 


0.98l°-! 7 7 


X— ray 


1 1 1 r>i i Ytr\^i"\r 
ILllllllH-JolLry 


12.60 


1.24 


0.72 


Abell2377 


0.081 


n oi+o:o5 

u - t5 — 0.05 


0-34l°° 6 6 


X— ray 


1 LlllllllUblLry 


12.07 


1.09 


0.28 


Abell2382 


0.062 


12+ 003 
— 0.03 


0.13l°° 3 3 


X— ray 


1 LlllllllUblLry 


11.43 


1.06 


0.30 


Abell2384 


0.094 


n ct;+0.06 

U.OO_Q Qg 


0.63l°-° 7 7 


X— ray 


ILllllllH-JiMLry 


12.59 


1.15 


0.23 


Abell2402 


0.081 


22+ ' 04 
u - zz _0.04 


24+ - 05 
u - z -0.05 


X— ray 


li i m i ti i~~\r 
i hi 11111 uoiLr^y 


11.67 


1.17 


0.41 


Abell2415 


0.058 


19+ 003 
— 0.03 


o.2iin 3 3 


X-ray 


luminosity 


11.46 


1.05 


0.26 


Abell2426 


0.098 


44+ ' 07 
—0 07 


o.5oin 8 8 


X-ray 


luminosity 


12.07 


1.02 


0.41 


Abell2428 


0.085 


0.251°'°* 


0.28in 5 5 


X— ray 


1 LlllllllU&lLry 


11.64 


1.07 


0.46 


Abell2443 


0.108 


O-SllHs 


0.39l°° 6 6 


X— ray 


1 LlllllllU&lLry 


11.98 


1.13 


0.53 


Abell2457 


0.059 


16+ ' 04 
u - ±u _0 04 


0.18l°;^ 


X— ray 


1 LlllllllU&lLry 


10.83 


1.03 


0.48 


Abell2495 


0.078 


0-29l° °4 


33+ 04 


X— rav 


luminosity 


11.69 


1.09 


0.37 


Abell2496 


0.123 


0.351°!! 


0-44l° : !! 


X-ray 


luminosity 


11.92 


1.17 


0.71 


Abell2589 


0.042 


0.20+°'°2 


0.22l°;° 2 


X-ray 


luminosity 


10.31 


1.04 


0.39 


Abell2593 


0.043 


0.14l°-°2 


16+ 002 

tJ *- LU — 0.02 


X-ray 


luminosity 


10.36 


1.02 


0.39 


Abell2597 


0.085 


0.38l°;° 7 


43+ - 08 
— 0.08 


X-ray 


^as mass 


12.31 


1.01 


0.25 


Abell2622 


0.062 


13+ - 02 
— 0.02 


u -- LO —0.02 


X-ray 


luminosity 


11.38 


1.01 


0.32 


Abell2626 


0.057 


21+ ' 02 
— 0.02 


n 9 o+0.02 
u -^°__ 0.02 


X-ray 


luminosity 


10.75 


1.09 


0.48 


Abell2627 


0.126 


n o 9 +0.06 
U ->J Z _0.06 


41+0-08 
u -0.08 


X-ray 


luminosity 


12.51 


1.22 


0.43 


Abell2717 


0.050 


12+ ' 02 
— 0.02 


14+ - 02 

— 0.02 


X-ray 


luminosity 


10.94 


1.06 


0.31 


Abell2734 


0.062 


26+ ' 03 
u - ZD -0.03 


29+ 03 
u - zy -0.03 


X-ray 


luminosity 


11.17 


1.05 


0.38 


Abell376 


0.049 


Uil -0.02 


0.l8l°;° 2 2 


X-ray 


luminosity 


10.78 


1.11 


0.35 


A hell 399 


0.072 


^2 +0 ' 05 
U - OZ -0.05 


0-59l°° 6 6 


X-ray 


luminosity 


10.84 


1.01 


0.70 


Abell401 


0.074 


I n+0.2 
1 -°-0.2 


1 5+°' 2 


X-ray 


^as mass 


10.91 


1.22 


0.69 


AbelU15 


0.197 


1 0+°' 2 


i q+0.2 


X-ray 


luminosity 


13.40 


1.25 


0.41 


AbelU201 


0.169 


n co+0.09 
U -°'J_0.09 


0.68t°;g 


X-ray 


luminosity 


13.16 


1.29 


0.40 


Abelll204 


0.171 


n KQ+o.09 


0.76t8;g 


X-ray 


luminosity 


13.58 


1.12 


0.28 


Abelll246 


0.190 


62 +0 1 ° 
u - DZ -o.io 


79+ - 13 


X-ray 


luminosity 


13.67 


1.32 


0.30 


Abelll423 


0.213 


-, 9 +0.3 
X - -0.3 


1 s+ - 4 


X-ray 


^as mass 


13.65 


1.41 


0.38 


Abelll553 


0.165 


5O+0.09 
U - OM _0.09 


76 +an 
u -' D -o.n 


X-ray luminosity 


12.60 


1.11 


0.64 


Abelll682 


0.234 


11 '-0.4 




X-ray 


^as mass 


13.12 


1.31 


0.72 
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Table Al. — continued from previous page 







Cluster 








BCG 




Name 


Redshift 


Mass 


Mass today 


Mass Proxy 


Ks 


J-Ks 


Stellar Mass 






hnl5 m l 


Tin 15 A//U1 

[1U 




[mag] 


[mag] 


[1U iVlQ\ 


Abelll704 


0.221 




92+ - 17 
u - az _0.17 


X-ray luminosity 


13.56 


1.23 


0.43 


Abclll758 


0.279 


o-87l° ;l 4 4 


1 q+0.2 


X-ray luminosity 


13.96 


1.38 


0.45 


Abelll763 


0.223 


9 q+0.5 
z -°-0.5 


q c+0.7 


X-ray gas mass 


13.11 


1.33 


0.67 


Abelll835 


0.253 


i 7 +0.2 
1 -'-0.2 


9 c+0.3 
z -°-0.3 


X-ray gas mass 


12.92 


1.44 


0.99 


Abelll914 


0.171 


1 4+0-2 


i-»±8:i 


X-ray gas mass 


12.85 


1.24 


0.54 


Abelll961 


0.232 


n tc+0.10 

u - 00 -o.io 


o-silaw 


X-ray luminosity 


13.52 


1.34 


0.49 


Abcll2009 


0.153 




Q1+ 013 
u - ai -0.13 


X-ray luminosity 


12.83 


1.19 


0.46 


Abcll209 


0.209 


-i 7+O.3 
i -'-0.3 


2 6 +0 ' 4 


X-ray gas mass 


13.07 


1.39 


0.62 


Abell2111 


0.229 


1 l+°' 2 




X-ray gas mass 


13.75 


1.35 


0.39 


Abell2163 


0.203 


5 2+°' 7 


8-211:1 


X-ray gas mass 


13.24 


1.74 


0.51 


Abcll2218 


0.176 


QfJ+016 
u - aD -0.16 


1 q+0.2 
i - J -0.2 


X-ray gas mass 


13.35 


1.11 


0.36 


Abell2219 


0.226 


n c+0.3 
Z -°-0.3 


3-9l2' 5 5 


X-ray gas mass 


13.34 


1.36 


0.55 


Abell2254 


0.178 


62+ - 09 
u - DZ -0.09 


0.80l° .l 2 


X— ray luminosity 


13.19 


1.27 


0.42 


Abell2261 


0.224 


1 - a -0.4 


o q+0.6 


X-ray gas mass 


12.62 


1.42 


1.06 


Abell2445 


0.165 


n o 7 +0. 08 
U -J'_0.08 


47+0- 10 
u -*'-o.io 


X-ray luminosity 


13.22 


1.18 


0.36 


Abell2561 


0.163 


n 09+0. 08 
u - oz -0.08 




X-ray luminosity 


13.53 


1.18 


0.27 


Abcll521 


0.248 


i-5l° : 2 2 


9 q+0.4 

z - J -0.4 


X-ray gas mass 


13.53 


1.33 


0.55 


Abell586 


0.171 


0.82l° ;l 2 2 


1 1+ ' 2 


X-ray luminosity 


13.13 


1.23 


0.42 


Abell661 


0.288 




1 c+0.3 
i -°-0.3 


X— ray luminosity 


13.53 


1.31 


0.71 


Abell665 


0.182 


1 7+0.2 
L -' -0.2 


9 9+O.3 
z - -0.3 


X-ray gas mass 


13.69 


1.23 


0.28 


Abcll68 


0.255 


1 n+°- 2 


1 c+0.3 


X-ray gas mass 


13.51 


1.43 


0.58 


Abell750 


0.180 


u -' z -0.12 




X-ray luminosity 


13.05 


1.31 


0.49 


Abcll773 


0.217 


1 2 +01 


1 7+0.2 
1 - ' -0.2 


X-ray gas mass 


13.22 


1.42 


0.58 


Abell907 


0.153 




0-8llo.il 


X-ray luminosity 


13.18 


1.32 


0.33 


Abell963 


0.206 


u - ai -0.13 


1-412-.2 


X-ray gas mass 


12.94 


1.37 


0.68 


RX J1720.1+2638 


0.164 


1 1+ 0,1 
1 - 1 -o.i 


1 4+01 
^ -0.1 


X-ray luminosity 


12.97 


1.21 


0.45 


RX J2129.6+0005 


0.235 


1 +0 ' 2 


1 k+0.3 
i -°-0.2 


X-ray gas mass 


13.27 


1.37 


0.63 


Zwl432 


0.186 


46+ ' 11 


59+0- 14 
u - oa -0.14 


X-ray luminosity 


13.31 


1.31 


0.41 


Zwl693 


0.225 


61+ - 15 


0-89l° :23 


X-ray luminosity 


13.24 


1.38 


0.60 


Zwl883 


0.194 


c 4 +0.13 


o-69±8;i; 


X-ray luminosity 


13.22 


1.28 


0.48 


Zw2089 


0.230 


U - 4Z -0.05 


o-60±°;8| 


X-ray gas mass 


14.10 


1.40 


0.28 


Zw2379 


0.205 


n 40+010 

u -^ a -0.10 


0.7ll° .l 5 5 


X-ray luminosity 


13.86 


1.27 


0.29 


Zw2701 


0.214 


54+0 09 
u -°^_0.09 


0.7812;! 4 


X-ray gas mass 


13.40 


1.30 


0.48 


Zw348 


0.255 


n 7K+0-16 
u - ' D _0.16 


l-llo.2 


X-ray luminosity 


13.78 


1.43 


0.45 


Zw3916 


0.206 


n c: 4 +0.08 
u - o4 -0.08 


o-78i2:ll 


X-ray luminosity 


13.80 


1.31 


0.31 


Zw5247 


0.195 


1 1+0.3 
^ -0.3 


1 4+0-3 


X-ray gas mass 


13.41 


1.28 


0.40 


Zw5768 


0.266 


86+ - 14 


1 0+0.2 

i - J -0.2 


X — rav 1 1 1 m 1 ni~it;i1"^7 


13.08 


1.25 


0.93 


Zw7215 


0.292 


0.84l°;l 7 


1 2+ - 3 
!- z -0.3 


X— ray luminosity 


14.12 


1.45 


0.43 


Abell2390 


0.228 


q 9+O.6 
-0.5 


5 0+ - 9 
°- U -0.8 


X-ray temperature 


13.49 ±0.07 




0.49 


MS0440+02 


0.197 


i-4l°;I 


1 8+ 1 - 1 
1 -8_0.4 


X-ray temperature 


13.34 ±0.05 




0.44 


MS0451+02 


0.201 


0.78l° ;f 5 


1 n+0.3 
L - z -0.2 


X-ray temperature 


13.94 ±0.07 




0.26 


MS0839+29 


0.193 


n 09+0.05 
u -° -0.05 


41+ 07 


X-ray temperature 


13.41 ± 0.06 




0.40 


MS1006+12 


0.261 


1 +0 ' 3 


1 5 +0.5 


X-ray temperature 


13.79 ±0.07 




0.47 


MS1231+15 


0.235 


46+ 14 


66+ - 21 

U.DD_ 21 


Velocity Dispersion 


13.89 ±0.07 




0.36 


MS1455+22 


0.257 


71+0- 25 

u -' 1 -0.16 


1 +0 ' 4 


X-ray temperature 


13.56 ±0.06 




0.57 
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Table A2. The intcrmediate-rcdshift subsamplc 







Cluster 








BCG 




Name 


Rcdshift 


Mass 


Mass today 


Mass Proxy 


Ks 


J-Ks 


Stellar Mass 






[1U 


[1U JHQJ 




Yin o rr 

[magj 


Yin o rr 

[magj 


[1U iVlQ\ 


MACS J0018. 5+1626 


0.541 


i 7 +0.5 
1 -'0.5 


, Q+1.3 

°- a 1.3 


X-ray temperature 


15.35 


1.52 


0.50 


MACS J0025.4-1222 


0.478 


95+ - 21 


4 -9 + 4 4 


X-ray temperature 


15.70 


1.67 


0.28 


MACS J0257.6-2209 


0.504 


2 1+ 0,4 


5.0++ 2 


X-ray temperature 


14.77 


1.61 


0.74 


MACS J0404. 6+1109 


0.358 




2 0+ 2 ' 


X-ray temperature 


13.82 




0.85 


MACS J0429.6-0253 


0.397 


i 4+O.6 


2 4+ 1 ' 1 


X-ray temperature 


13.58 




1.34 


MACS J0454. 1-0300 


0.550 


1 0+ ' 3 
1 - u 0.3 


9 O+0.8 
z -°0.8 


X-ray temperature 


15.29 


1.58 


0.54 


MACS J0647.7+7015 


0.584 


o k+0.5 
z -°0.5 


6-0+i' 2 


X-ray temperature 


14.87 


1.76 


0.89 


MACS J0744. 8+3927 


0.686 


1 1+ ' 2 
1 - ± 0.2 


3.0+° 5 8 


X-ray temperature 


15.33 


1.88 


0.76 


MACS J2129.4-0741 


0.570 


i q+0.2 
1 ' Z 0.2 


9 7+O.6 
Z -'0.6 


X-ray temperature 


15.57 


1.72 


0.45 


MAPS 7994^ 0+9fi37 


0.301 


60+ - 14 






14.03 




0.49 


MS0016+16 


0.547 


1 6+°' 4 


q 7+1.0 
°-'-0.8 


X— ray temperature 


15.29 + 0.08 




0.54 


MS0302+16 


0.425 


n 9Q+0-65 
U - ZM -0.13 


0.53+1-30 

— U.25 


X-ray temperature 


15.01 + 0.07 




0.42 


MS0451-03 


0.539 


2 0+ ' 4 


4 6+ 11 


X-ray temperature 


15.18 + 0.07 




0.58 


MS1008-12 


0.306 


70+ ' 24 
u - ' u -0.16 


1 2+ ' 4 
1 - Z -0.3 


X-ray temperature 


13.68 + 0.08 




0.70 


MS1224+20 


0.326 


75+0- 25 


i q+0.4 


Velocity Dispersion 


14.41 ± 0.08 




0.41 


MS1358+62 


0.329 


i 7+0.9 
1 -'-0.4 


2 9+ 1 ' 6 


X-ray temperature 


14.29 + 0.06 




0.46 


MS1512+36 


0.373 


21+ 012 


34+0 20 
U- 13 -0.13 


X-ray temperature 


14.63 ± 0.08 




0.44 


MS1621+26 


0.427 


94+104 


1 S+ 2 ' 2 
1 - 8 -1.0 


X-ray temperature 


14.98 + 0.07 




0.44 
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Table A3. The high-rcdshift subsamplc 







Cluster 










BCG 




Name 


Redshift 


Mass 


Mass today 


Mass Proxy 


Ks 


J-Ks 


Stellar Mass 






fin 15 Mr-A 

[1U 


[1U 




[magj 


[mag] 


[1U 


SpARCS J003442-430752 


0.867 


n o fi +0.16 


i 1+0.6 
'-■'--0.6 


Velocity Dispersion 


16.52 ±0.04 


1 


.86 ± 0.03 


0.38 


SpARCS J003645-441050 


0.867 


n 45+016 




Velocity Dispersion 


16.09 ±0.05 


1 


.84 ± 0.03 


0.56 


SpARCS J161314+564930 


0.873 


2 6+0' 6 




Velocity Dispersion 


15.69 ±0.01 


1 


.79 ±0.01 


0.81 


SpARCS J104737+574137 


0.956 


u - zy -0.13 


i oo+o- 40 

1.UU_q 4 g 


Velocity Dispersion 


17.14 ±0.03 


1 


.89 ±0.03 


0.24 


SpARCS J02 1524-034331 


1.004 


26+0- 17 
U - ZD -0.13 


1 0+0' 8 
^-o.s 


Velocity Dispersion 


16.88 ±0.14 


1 


.86 ± 0.05 


0.33 


SpARCS J105111+581803 


1.035 




n 42+ 012 

U - 4Z -0.23 


Velocity Dispersion 


16.88 ±0.05 


1 


.74 ± 0.03 


0.35 


SpARCS J161641+554513 


1.156 


0.281°;^ 


i n+0.6 

1 '°-0.6 


Velocity Dispersion 


17.02 ±0.03 


1 


.73 ± 0.02 


0.36 


SpARCS J163435+402151 


1.177 


n 44+0- 11 


9 i +0.6 
-0.9 


Velocity Dispersion 


17.35 ±0.02 


1 


.84 ± 0.03 


0.27 


SpARCS J163852+403843 


1.196 


n 10+ 003 


0.38±°;£ 


Velocity Dispersion 


17.65 ±0.05 


1 


.91 ± 0.05 


0.21 


SpARCS J003550-431224 


1.340 


oq+0.13 


2.5+.J;? 


Velocity Dispersion 


17.52 ±0.01 


1 


.98 ± 0.01 


0.29 


SpARCS J033056-284300 


1.620 


94+0' 10 




Velocity Dispersion 


17.88 ±0.04 






0.29 


SpARCS J022426-032331 


1.630 


04+ 01 


26+ 05 
U - ZD -0.17 


Velocity Dispersion 


18.07 ± 0.03 






0.25 


CL J0152. 7-1357 


0.830 


n 45+0- 17 


1.4+0-6 
— 0.5 


X-ray temperature 


16.96 ±0.08 


1 


.80 ±0.08 


0.23 


XLSS J022303. 0-043622 


1.220 


n , .+0.04 


0.68^0 


X-ray temperature 


17.72 ±0.01 


1 


.82 ± 0.01 


0.21 


XLSS J022400. 5-032526 


0.810 


n 19+005 




X-ray temperature 


16.49 ±0.10 






0.34 


RCS J0439-2904 


0.950 


o.o3ir 01 


0.08±°;°* 


X— ray temperature 


17.70 ±0.08 


1 


.86 ±0.08 


0.14 


2XMM J083026+524133 


0.990 


0.981^ 


4-otl;J 


X-ray temperature 


16.58 ±0.05 


1 


.90 ± 0.06 


0.43 


RX J0848.9+4452 


1.260 


n 47+0- 17 
u -*'-0.l3 


2 7+ 1 ' 2 


X-ray temperature 


17.00 ±0.02 


1 


.86 


0.42 


RDCS J0910+5422 


1.110 


n tre+0.30 
U - OO -0.19 


9 7+1.8 


X-ray temperature 


17.88 ±0.05 


1 


.83 ± 0.06 


0.15 


CL J1008.7+5342 


0.870 


n 19+010 


n 54+o- 32 

u -°^-0.19 


X-ray temperature 


16.42 ± 0.08 


1 


,97 ± 0.09 


0.41 


RX J1053.7+5735 West 


1.140 


n nc+0.04 


1 l+°' 2 


X-ray temperature 


17.21 ± 0.06 


1 


.99 ±0.07 


0.29 


MS1054.4-0321 


0.820 


n 97+0 28 
u - M '-0.22 




X— ray temperature 


16.04 ±0.10 


1 


,80 ±0.10 


0.53 


CL J1226+3332 


0.890 


1 s+0- 4 




X-ray temperature 


16.00 ±0.06 


1 


,71 ± 0.07 


0.63 


RDCS 11959 9 9997 


1.240 


66+ 08 


o q+0.6 


j\- lay Lt_.ui ljui ai u.1 


J- 1 .ou in u.uo 


1 


,83 ±0.01 


0.30 


RDCS J1317+2911 


0.810 


24+0 .19 
u - z -0.09 


71 +0.66 
u - ' -0.29 


X-ray temperature 


17.27 ±0.15 


1 


,68 ±0.17 


0.17 


WARPS J1415. 1+3612 


1.030 


54+0- 15 


9 n+0.8 

A - 6 -()Ai 


X-ray temperature 


16.76 ± 0.04 


1 


.86 ±0.05 


0.39 


CL J1429.0+4241 


0.920 


57+0- 32 
U - ' -0.17 


2 1+ 1 ' 4 


X-ray temperature 


17.43 ±0.20 


1 


,78 ±0.22 


0.18 


CL J1559.1+6353 


0.850 


25+ ' 21 

u - zo -o.ii 


74+0- 72 
U - ' 4 -0.35 


X-ray temperature 


17.21 ± 0.09 


1 


,90 ±0.12 


0.19 


CL1604+4304 


0.900 


09+0- 10 


n 9 o+0.30 


X— ray temperature 


17.61 ± 0.09 


1 


,68 ±0.12 


0.15 


RCS J162009+2929.4 


0.870 


o si+o- 37 


n q5+ 133 


X-ray temperature 


17.63 ±0.12 






0.14 


XMMXCS J2215.9-1738 


1.460 


0.l8t°;°? 


1 ' 1 -0.5 


X-ray temperature 


18.72 ±0.01 


1 


,84 ± 0.02 


0.11 


XMMU J2235.3-2557 


1.390 


0.88±°;2° 


6 5+ 2 ' 7 
— 2.0 


X-ray temperature 


17.34 ±0.01 


1 


,87 ± 0.02 


0.37 



Table A4. BCGs without cluster mass measurements 



Name Redshift Ks J-Ks Stellar Mass 

[mag] [mag] [10 12 M ] 



Abell2292 





,119 


12, 


06 


1.06 





,59 


Abell2665 





,056 


10. 


74 


1.11 





,47 


Abell291 





,196 


14, 


10 


1.32 





.22 


MACSJ0150.3-1005 





.363 


13. 


.90 







.82 


MACSJ0329.6-0211 





.451 


14. 


13 




1 


,07 


MACSJ1359.8+6231 





,330 


14. 


32 







.45 


MACSJ2050.7+0123 





,333 


14. 


67 







.33 


MACSJ2214.9-1359 





.495 


14, 


71 


1.67 





.76 


MACSJ2241.8±1732 





,317 


14, 


39 







,39 


RCS0224-0002 





,770 


16, 


87 







,22 
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